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AERODYNAMIC MATHEMATICAL MODEL 


Definition of the Problem 


To study active control of the aeroelastic response of rotor blades 
to gust or discrete pulse loadings, an analytical aerodynamic model is 
needed which allows concentrated loadings to be applied in a general way. 

The model should also be in a form suitable for use with the electro- 
mechanical analog used in the experimental simulation study. 

Aerodynamic strip theory is a convenient way to formulate a model 
that allows for general spanwise and timewise loads. To use strip theory, 
three-dimensional effects (for example, spanwise aerodynamic derivatives) 
must be neglected. The basic relations between the unsteady blade motions 
and the unsteady air flow are derived from quasi- steady aerodynamic theory. 
This is equivalent to assuming that the actual unsteady, Theodor sen-like 
aerodynamic load functions are replaced by unity in the lift and moment 
equations of atrip theory (Ref. 1). There is a controversy in the literature 
about the aerodynamic damping coefficient derived by this approximation 
(Ref. 2), but the approximation is certainly used widely. The unsteady 
effects of the rotor wake are ignored, which is a necessity in any case be- 
cause a tractable description of the wake aerodynamics is lacking. 


Specific Assumptions 

In addition to the use of quasi- steady strip theory, a number of 
specific assumptions are made in developing the mathematical model. For 
convenient reference, these are listed below as a group. 

1. Spanwise velocities do not make a direct contribution to 
lift or drag, 

2, The induced downwash and the component of the forward 
flight velocity normal to the rotor disk are both small rela- 
tive either to the tip rotational speed or to the forward 
velocity. 

3. The shaft tilt angle is small, so that sin ttg = ttg and 

cos ttg = 1. 

4, The bending, flapping, and torsional deflections are so small 
that products of them with other small quantities are negligible. 
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5, Cyclic pitch 0^ is small, so that sin = 0(, and 
cos 9c = 1. 

6. The slope of the lift-coefficient versus Uj. curve for re- 
versed flow is the same as for forward flight, for the 
unstalled flight regimes. 

7, Blade stall is confined to blade sections and azimuth angles 
where the chordwise velocity is close to zero. Loads due 
to stall effects, therefore, can be neglected because of the 
smallness of the dynamic pressure in the stalled region. 

8. The rate of climb is zero. 

The aerod3mamic assumptions are justified in Reference 3. 
Assumption No. 6 is a good approximation for typical symmetric airfoil 
sections (such as NACA 0012) used for helicopters, when the angle of 
attack is below stall. With regard to Assumption No. 7, the angle of 
attack at a section varies from non- stalled positive values to zero, and 
then very rapidly through the range of large values of angle of attack 
(i. e. , stall) as the chordwise velocity varies from a small forward value 
to a small rearward value; in the region where the chordwise velocity 
has a large rearward value, the angle of attack is near -tt and the blade 
section is again unstalled according to Assumption No. 6. 

Ih the light of these assumptions, no pretense is made that the 
model is "exact" or likely to be of value for flutter analyses and large- 
amplitude or large angle of attack aeroelastic instability studies. The 
model is suitable, however, for a general study of active control of rotor 
blade response by systematic variation of the parameters, the purpose 
for which it is formulated. 


Model Development 

Basic Equations : Figure 1 shows schematically the blade under 
consideration. All the velocities are depicted relative to the blade. The 
mathematical development is somewhat similar to that presented in Refer- 
ences 3, 4, and 5, but without the restrictive assumptions inherent in 
Reference 5. Since the present study includes elastic motions of the blade 
as well as flapping and feathering, there is no particular advantage to any 
choice of reference axes for the angular quantities; therefore, the axis of 
the rotor shaft is selected herein. With regard to standard terminology 
(Reference 6), the shaft angle is tilted from the vertical through the angle 
ttg with positive ttg conventionally being taken as rearward. The rest 
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o£ the terminology corresponds as much as is possible with Reference 6, 
except that the blade- section pitch angle 9 is the angle between the line 
of aero lift of the blade section and the plane perpendicular to the shaft 
axis, rather than the plane perpendicular to the axis of no feathering. 

The chordwise velocity at the quarter- chord point, in the plane 
perpendicular to the rotor axis is 

U<p = fix + V sin fit + fw - (Zga-ZjjQ) 0 + (2^) Y 

- (Za) (3w/'dx) 0 cos 0 } sin 0 (1) 

where cos ttg has been replaced by unity in the second term on the right. 
The term (Zi) fdw/3x) cos 9 is derived with reference to Figure 2a. 
Because of tbe downward deflection of the blade, the rotation vector 
(oriented parallel to the Yq axis) has a component parallel to the elastic 
axis (that is, in the plane of the blade) of magnitude fi sin[tan“^(3w cos ®/3x)]. 
Since cos $ = cos(0 + y)> neglecting products of small terms gives 
n(3w/Bx) cos 9 for this rotation component. The velocity in the plane 
perpendicular to the rotor axis, due to this rotation about the elastic axis 
is Zi times this angular velocity. Terms such as this have been variously 
called gyroscopic or Coriolis effects in the literature. 

The velocity at the quarter- chord point parallel to the rotor axis and 
directed upward is 

Up = Vttg - V + (w - (Zea-Z^o^ 0+ (Zj) Y 

-- (ZJ) (Bw/Sx) Q cos 9} cos 9 V (5w/Bx) cos 9 cos fit (2) 

where sinag has been replaced by in the first term on the right. 

The origin of the term (Zi) (3w/3x) n cos 9 has been explained previously. 
The last term on the right Is derived with reference to Figure 2b, Be*- 
cause of the blade deflection, the spariwise velocity V cos Ctg cos fit can 
be resolved into components parallel to and perpendicular to the blade sur- 
face. The component parallel to the span can be neglected {Assumption 
No, 1) and the pei^pendicular component can be thought of as having a com- 
ponent along the Yq axis equal to the spanwise velocity times 
sin [tan"^(Bw/5x) cos®] times cos[(&w/hx) cosi], which is equal to the 
term in question after the small angle assumption is made. 

The total velocity U or ^^relative wind^^ at a blade cross-section- 

r 2 21 i 

is LUt + Up J jr and the angle made by U and U«j' is the upwash or 
inflow angle 0, as shown in Figure 1* The true angle of attack at a 
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cross-section is = 5 -f 0, which is the angle between the chordline 
of a syimnetrical foil and the relative wind. 

Care must be exercised in expressing the lift and drag coefficients 
for the blade, because of the possibility of reversed flow for x < V/fi. 

In developing the lift coefficient relations, only the unstalled regions of 
angle of attack for forward flight (around ttj, = 0) or for reverse flow 
(around CL^= -rr) need to be considered, because of Assiimption No. 7. 

For unstalled forward flight, the lift coefficient is = Clq ^'■r "where 

Clo slope of the lift coefficient versus o.^ curve near ttj. = 0. 

For unstalled reverse flow, the angle of attack is close to -tt, and the 
lift coefficient is -0^0 ■** positive lift direction is still upward, 

which does not coincide with the conventional definition for which positive 
lift is obtained by rotating the relative wind through tt/2 clockwise. 

(This slight departure from conventional practice is desirable because 
it keeps the direction of positive lift upward over the entire span. ) When 
U>p > 0, and Up is small (forward, unstalled flight), the inflow angle 
0 s tan“^ Up/U>p is well approximated by 0 = Up/U-p ; when U>p < 0 
and Up is small (reverse, unstalled flight), the inflow angle is well 
approximated by Up/U>p + n, i. e. , the inflow angle is near n when 
referred to the positive Zq direction. Thus, the lift coefficient for a 
blade section in unstalled forward flight is Cp, = Clq (0 + Y + Up/UT) 
Sind for unstalled reverse flight it is -Cj^q (0 + Y + Up/Up) . The lift 
arising from unsteady pitching, flapping, and bending for the quasi-steady 
approximation is given by (Reference 1) 

Fi = ttp(c/2)2 {w - (Zea-Zxo) ^ + (^i) Y 

+ U[0+ Y- Ow/3x) cos 9]} (3) 


The small "added-mass" terms can be neglected, so 
Fj = ^TTpc^ U[ 0+ Y- (3x/Bx) cos 0] 


(4) 


Combining Equation (4) with the angle -of-attack contributions gives the 
total lift as 


Fl = 4 pc U^ Clo 


_Turf, 


(e+ Y+ up/0x) 


+ ^npc^ U [9 + Y“ 0 (Bx/Bx) cos 9] 


(5) 


The term j^UT/IUTl changes the sign of the angle -of- attack lift contribu- 
tion in the reversed flow regime, as is required. By using the expressions 
for U-p and Up given in Equations (1) and (2), and neglecting products 
of small terms, the section lift is derived as 
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j(6+ y) 

- pc CjjQ (Clx+V sin fit) I^V- Vctg - w cos 0 - ® 

+ xnpc^ {flx + V sin fit) 0 + Y “ fi cos bJ (6) 


To simplify the analysis, it has been assumed in Equation {6) that the 
elastic axis, the cyclic pitch axis, and the quartei*- chord point coincide 
(Zga = ZjjQ, Zi = 0); this assumption is not a necessity but does reduce 
the number of parameters that need to be varied. 


The drag force is derived similarly; it is 


■r, 1 ^ ,r - fix+Vsinfit 

Pd = *p= Coo(fi=:+ 


(0 + Y) 


- pc Cj 3 Q (fix + V sin fit) j^v - Vttg - w cos 9 


(7) 


Here, CpjQ is a computed slope derived from the drag coefficient versus 
curve, such that the drag coefficient has approximately the correct 
values over the operating range of tty (say, 4*^ to 12°). This procedure 
cannot be extended to very small values of tty because of the non- zero 
value of drag coefficient for tty = 0; however, it is almost a necessity 
that the drag force and the lift force have similar functional forms for use 
with the electromechanical analog, in order to avoid undue complexity. 

It is also assumed that Cqq is the same for the unstalled forward and 
reverse flow regimes. This latter is a fairly crude assumption, but the 
drag force is itself small for unstalled flight, so the absolute error in 
computing the thrust is tolerable. The total section thrust T is equal 
to Fjj cos - Fjj sinay. Thus, the final expression for the thrust per 
unit span at a blade section after neglecting products of small terms is 


T = ^ pc Cijq(co3 9- Cj 5 q/ Clq sin 0) (fix+ V sinfit)^ X 
fix+Vsinfit1 _ , „ _ 

jfix+VsinfitiJ Y) - s pc Clo(co8 9- Cdq/Clo 9) 


(fix+ Vsinfit) 


V - Veto - w cos 0 - V 



cos fit cos 9 


h^TTpe^ cos 0(fix+ Vsin fit) ^0 + y - 0 (^) 




(8) 
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The aerodynamic torque about the rotation axis caused by the lift 
and drag on the blade section is 

□x = X (Fl sicMS + cos « (Fl 0 + Fd)x (9) 


Thus, after neglecting products of small terms, the section torque turns 
out to be 


Qx = i P<^ *^110 V sin Qt) 

- V cos Qt cos bJ 6 

+ ^ pcxCj^Q (Qx+ V sinQt)^ 


V - Vap! - w cos 9 


Ox+ V sin Qt 
fix+ VsinCltj 


(0+Y) 


( 10 ) 


Integrated Eqxiations ; In order to permit the electromechanical 
analog to be ”de-bugged, '* and to allow values of ttg, 0^,, and v to be 
selected with a minimum of trial-and-error, as well as *:o provide a quick 
means of verifying that the operational parameters selected are in the 
range of helicopter practice, the thrust and torque equations are integrated 
(approximately, for forward flight) over the blade span. The blade is 
assumed to have a built-in, linearly varying pitch along its span; 

0 = 0Q+ Qq (x- 6q), where 9g is the pitch at the blade root, 0j is the 
difference in the pitch between the tip and the root, and 9 q^ = Gj/X is the 
rate of twist. The cyclic pitch is assumed to have a constant amplitude 
along the entire span, i. e. , there is no dynamic, elastic twisting of the 
blade. 


For hovering conditions, the cyclic pitch, 0^, and V are both 
zero; for hover, also, the flapping and bending motions do not affect the 
time-averaged thrust, so Equation (8) can be integrated exactly to give 
the thrust per blade; 

Tavg = pc Clo o- 0o eo) ( 1 - E.^) + I R( 1 - Eg^) } 

-ipcCLO V (1 -Eq2) (11) 


where R= & + eg, Eg= Sg/R, cos 9 has been approximated by unity, 
and the small drag contribution has been neglected. For hover, simple 
momentum theory (Reference 6) is used to predict the induced velocity; 



(12) 


where b is the number of blades. The simultaneous solution of Equations 
(11) and (12) permit v to be computed for any desired value of thrust. 

The total torque for the hovering case is, from Equation (10), 
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( 13 ) 


^:otal= |(®o " ®o) j^sfiL ^ 4C]joj 

■ 8 / R f ^ . . .^JDO 

° [4fiL 5 CloJ 


(Equation (13) is used px-imarily to determine whether the predicted torque 
has a realistic value for the desired thrust. ) 


For the forward flight, the blade- section equations camaot be inte- 
grated exactly because many of the terms depend on the elastic response 
of the blade. However, approximate expressions can be formed and used 
for guidance. These are now developed. The induced velocity and the 
total thrust, averaged over one blade revolution, are related by the over- 
all momentum equation (Reference 3); 


v"^ + 2V v%g + V^V^. - 


'b T 


BIK 
L 2rrp . 


= 0 


(14) 


In this equation, the angle of attack of the rotor disk (tip-plane path) has 
been replaced by the shaft tilt angle ttg ; these angles differ by the cyclic 
pitch, but since only a first approximation is desired as a starting point, 
this simplification is acceptable. In fact, in most cases both the first and 
second terms in Equation (14) are negligibly small. Letting Vq be the 
solution when the first two terms are neglected, then 


V 


o 


^ "^avg 

Lettp r2. 


(15) 


The next approximation to the solution of Equation (14), for small values 
of ttg, can be derived by an expansion in powers of ttg , The result is 


V = Vq 


“ 3 (Vp/V) ~ 

. 1 + 2(Vo/V)2 _ 




( 16 ) 


Equation (16) is one relation between v and ttg, which are the ’'unknowns*' 
for a given thrust and forward velocity. Actually (9v/Sag) is always very 
small, so it is a good approximation to assume v = v^. Thus, the approxi- 
mations in Equ' iion (14) with respect to the rotor angle of attack is not a 
serious limitation. The true value of the shaft tilt angle ttg must be se- 
lected so that Equation (8) gives the correct value for the desired average 
thrust, when the cyclic pitch 9^. has been simultaneously selected to 
diminish to y-ero as nearly as possible the one- per -revolution flapping 
(elastic and rigid body components) of the blade tip with respect to the 


shaft axis. la px*actice, the one-per-revolution thrust component along 
the shaft axis would probably be minimiaed. However, both of these pro- 
cedures result in almost the same cyclic pitch setting, and minimization 
of the one-per- revolution flapping is more convenient for the electro- 
mechanical analog, in order to reduce the amount of trial-and- error in 
choosing the correct V - -ag values with the analog, the one-per- 
revolution variation in the thrust from Equation (8) is set equal to zero; 
this, in combination with Equation (15) or (16), gives a set of first 
approximations. The blade flapping (which should be nearly zero), the 
elastic motions, and the reversed flow region are neglected in deriving 
these relations, although, of course, they are retained in the analog 
development. Integrating Equation (8) along the span, with w = 0, neg- 
lecting reversed flow, and setting the one -per -revolution time-varying 
thrust terms equal to zero, the amplitude of the cyclic pitch can be 
deduced to be equal to 



- 3(9o-9o^eo) -9jR 
1 + 9/4 (i^ 


(17) 


ha this equation, |u = V/QR is the advance ratio. The phase angle, with 
respect to the downwind position, for the cyclic pitch (when the cyclic pitch 
is expressed as 0gcos(nt+ i&q) ) is also derived in this way: 


to 


-tan 


-1 


6ttc 


[SRClo (1 + 


f 


(18) 


Equating the time-average thrust for forward flight to the thrust for hover 
gives the needed relation between CL^, v, and V: 


V 

“s = y 


V 


+ U(9o-0jeo+t e^^R) (l+9/4|i2) -4M[3(0o-9oeo)+6o 

^ _ (19) 


1 + 3/4 


In this equation, Vfj is the induced velocity for hover with the same 
thrust. 


Equation (15) is used to compute the value of v used by the analog 
for forward flight. Equations (17), (18), and (19) are merely used as 
•'first guesses" for 0^,, i|)o, and Ug; the true values are then determined 

by trial-and-error on the analog until the desired thrust is obtained and the 
one -per -revolution variation in the blade tip flapping is miniraized. 
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MODSLr^TG CONCEPTS 


Genei'al 


Modeling of the dynamic response of an artictilated helicopter rotor 
blade in forward flight can be accomplished by means of a relatively simple 
dimensionless analysis which results in the following nondimens ional 
equation, 


^ ‘ t ( PcA £ fo JL JiL 9 F \ 

& " m ' ' nR' OR* ' m&^QZj 


( 20 ) 


•where all the dimensional variables have been identified in the previous 
section, except for an arbitrary force F. As previously mentioned, the 
electromechanical apparatus -with -which the simulation was to be accom- 
plished was already available foi* the most part. (The apparatus will be 
described in more detail in the next section. ) Therefore, in order to 
accomplish a reasonable simulation of a representative prototype aircraft, 
care had to be exercised in developing the corresponding protot 3 rpe and 
model variables. A brief description of this process will now be given. 


From Equation (20), it is obvious that d'ynamic responses for only 
flapping and bending modes of the rotor blade will be considered. A fre- 
quency plot of these modes for the existing model apparatus is shown in 
Figure 3. Measured frequencies are compared with values predicted from 
simple beam theory. Note that results are presented up to a maximum 
frequency to include the third bending mode (13 per rev) only. Above this 
range, the results become suspect because of force increment resolution. 
E-vidence of this will be discussed later. 


Thus, all geometric and structural parameters for the blade were 
already fixed, and simulation had to be based on a proper selection of 
blade rotational veloci'ty and air density. By trial and error, a reason- 
able choice for these parameters was made, and the res'ults are shown in 
Table 1(a). By the substitution of these results and fixed parameters of 
the model appara-tus into the nondimensional expressions of Equation (20), 
the corresponding protot-ype parameters were developed. They are given 
in Table 1(b), and are considered to be t-ypical for a representative heli- 
copter. A further demonstration of the appropriateness of these protot-ype 
values will be given in a later section. 

Vortex Impingement 

The conditions existing when a vortex shed by a preceding rotor 
blade impinges on the following blade are described by Scheiman and Ludi 
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(Reference 7) and Ward (Reference 8). Impingement occurs at the inter- 
section of the following blade with the path of the preceding blade tip. The 
trailing vortex of the preceding blade may significantly affect the vibra- 
tory blade loading of the following blade, producing harmonics of all 
orders, including significant effects at higher orders. 

To simulate the trailing vortex type of impulsive loading phenomena 
in the model, a time dependent signal, v(t), was summed in with the in- 
duced downwash velocity analog signal in the analog circuit portion of the 
electromechanical apparatus. This time- dependent signal was a square 
wave whose duration was initially chosen to simulate prototype vortices 
of 0.6 meters (2.0 ft) in diameter and 8,5 m/sec (27,8 ft/aec) in ampli- 
tude. Later, other values of duration and amplitude were used (see 
Table 4). 

Three types of conditions were simulated in terms of the occur- 
rence of vortex impingement on the following blade. The first condition 
was a single pulse applied only once, and the second was a repeated 
"one -pel* -rev" condition in which v(t) was a series of pulses spaced at 
time intervals equal to one revolution of the blade. The pulses were 
introduced to the analog circuit controlling the exciter location at the 
point 88% of the distance from the center of rotation of the blade to the 
blade tip. The third condition was a "moving pulse" condition in which 
v(t) was a series of pulses, with a given time interval separation, all 
occurring within one blade revolution. Each pulse in the series was intro- 
duced to the analog circuit controlling different exciter locations along 
the blade (see Table 4 for locations). The series of pahves was repeated 
every revolution of the blade. 

The time interval separation of the series of pulses for the moving 
pulse condition simulated the situation in the prototype where various 
points along the following blade intersect the path of the preceding blade 
tip at different times. The time of vortex impingement for any location 
on the blade depends on the forward flight speed of the aircraft, or advance 
ratio, |i, the number of blades, b, and the rotational velocity of the 
blades, fl. 

The relevant equation is: 

tan [fltg - ?] = tan Qtj - (V/x) (t2 -tj) (21) 

where 5= 2rr/b and b is the number of blades, V is the forward flight 
velocity, x is the location on the blade, and tg and tj are the elapsed 
times of travel of the following and preceding blades, respectively. 
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Figure 4 {obtained by solving the above equation) shows how the 
point o£ blade interception by the vortex varies with time for two values 
of forward flight advance ratio. The ji = 0, Z6 type movement was not 
incorporated into the analog circuits because of complexity of implemen- 
tation, Thus, only monotonic type variation in impulse loading movement 
was studied. For the li = 0.42 prototype case, the time interval be- 
tween pulses for the monotonically moving pulse condition was obtained 
by locating the corresponding abscissa (time of vortex impingement) for 
the given ordinates (locations on the blade, i. e. , exciter locations) from 
Figure 4 for the 0.42 advance ratio. This is a series of pulses which 
moves with time down the blade from the blade tip towards the root, hence, 
a "moving pulse. " 

The tip vortex type simulation with monotonic movement is only 
valid for the prototype, 4-bladed rotor case, at an advanced ratio of 0.42. 
However, the monotonic variation was also used for the hover case and the 
|j = 0. 26 case. This provided additional dynamic response data that re- 
flects a decrease in number of prototype blades, in the case of |j = 0,26, 
and also provides an indication of blade response to an impulse, such as 
an imposed shock front or straight line vortex, moving across the rotor 
disk in hovering flight. 


ELECTROMECHANICAL ANALOG SIMULATION 


System Description 

Modeling of a lielicopter rotor blade experiencing the in-flight en- 
vironment was accomplished in this stndy by use of an electromechanical 
apparatus similar to that described in previous work {References 9 and 10). 
This apparatus consisted of a pin- free blade rotating at a constant angular 
velocity and a group of specially- developed analog circuits. Details of 
the mechanical apparatus, including the blade and rotational device, are 
described in References 9 and 10. The analog circuit is described below. 

The salient features of the blade geometry may be seen in Figure 5. 
The blade is pinned at the flap hinge a distance eo from the center of 
rotation. The response of the blade at various points (V 2 through VIO) 
could be measured from the output of velocity- sensing coils located at 
those positions. The blade was excited by force coils at the exciter posi- 
tions shown in the figure. The size and complexity of the analog circuit 
required to generate current for the force coils limited the number of 
forcing points which could be used to excite the blade. At many locations 
two adjacent force coils were coupled in parallel to extend the length over 
the blade to which the force from one circuit was applied. Where coils 
were paralleled (FI, F3, F4, F5) the force was assumed to be applied 
over a distance of 19. 1 % of the length from the tip of the blade to the 
center of rotation and aro\ind the center point between the exciter loca- 
tions, For the single coil, F2, the force was applied around point F2, 

9 . 6 % of the tip to rotor axis distance. 

Generation of the current to drive the force coils at the exciter loca- 
tions was accomplished with an analog circuit and associated apparatus. 

As an example. Figure 6 shows a block diagram of the part of this appar- 
atus used to generate the force at location F5 (similar circuits “/ere used 
to generate the forces at the other excitation points). Input to the analog 
circuit was the blade response at V 8 and VIO. These responses were 
summed and divided to obtain a velocity signal equal to their amplitude 
average at any time, W 5 . This velocity signal was assumed to be equiva- 
lent to the velocity response of the beam at midpoint for the excitation 
points of F5, X 5 .* The difference between the blade response signal from 
VIO and V 8 was integrated with time to obtain a signal propostional to 
the instantaneous difference in displacement of the blade, W5, at locations 
VIO and V 8 . Thus, a signal Ow/3x)g, the displacement slope of the 


■^Note that the V9 signal could also have been used for this response; how- 
ever, the use of V 8 and VIO was required anyway in developing the slope 
of displacement in the forcing increment. 
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blade at could be generated. Constants input to tbe analog circuit 

•were: d. c. voltages simulating the time invariant value of blade pitch. 


0 g , at blade increment 5; the induced downwash velocity, 
product of the for'ward flight speed and the shaft tilt angle. 


V 


o » 


va 


s ' 


the 

and 


the product of the rotation speed and xg, fixg. A one-per-rev generator 


attached to the mechanical simulator rotor axis generated a signal, 
sinOt, input to the analog circuit, and a phase shifter was used to ob- 
tain sin (fit + Iq) as another input. 


The details of the many analog computations are shown only in 
rough form in the block diagram. The last element of the circuit was a 
power amplifier (voltage-to- current converter) used to convert the signal, 
F5 into a current to drive the force coils. Referring back to Equations 
( 6 ) and ( 7 ), the equation for F5 the oxitput of the analog circuit, in terms 
of circui-t input signals is;"^ 


F5 = -0.25tt Axgpc^ cos 0g (flxg + V sin fit) 


“ (i) 


cos 9g 


[ JJU 

- 0.5 Axgpc Cj_|Q fcos 0g - 


sin 65 ^ (fiX 5 + V sin flt)^ 


Un 


IutI 


+ 0.5 Axgpc Cj^Q [cos 9g 


C 


DO 


C 


LO 


sinGg^ (nxg + VsinQt) X 


Vq - Vttg + v(t) - -Wg cos 65 - v 


/^\ 

1,9xA 


cos fit cos 9 




( 22 ) 


where 


Axg 

= 19 . 1% (eo t Ji) 

^5 

00 

00 

(e,3 + i) 

•Wg 

= (VIO 

+ V8)/2 

(9w/3x)g 

= J(V10 

- VS)/Axg 

CD 

= Gq + 

®o (^ " ®o) 

65 (t) 

= 0g + 

Gj. sin (Qt + 


e. 


= 0Q cos (fit+ IITq) 
and 


and where 9^ and are cyclic pitch amplitude and phase, respectively. 

0Q^ is the blade twist per unit length; Jt, is the blade length; p is the air 
density; c is the blade chord; and Cjjo ^LiO slopes of the 

drag and lift coefficients versus angle of a'ttack for a blade shape NACA 0012. 


‘Note that in Equation (22), the algebraic sign has been set to minus the 
aerod' 3 mamic force. This is done in the analog in order to make a positive 
force correspond to a positive w, to conform with ordinary beam theory. 
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Values for tlie fixed parameters for this study: JJ, c, eo> P» 0]^, 

b, CjjQ, and Cjjo given in Table 1. Values for the parameters 
varied during the study: Sq* V, ttg, ilfo* given in 

Table 2. The parameter v(t) was generated by a pulse generator as 
described in the section on Vortex Impingement. 

Procedures for initial Conditions 

The electromechanical model was run for a series of tests to simu- 
late in-flight conditions of a helicopter rotor blade. Those conditions for 
various advance ratios, |a, in which v{t) was zero (no vortex impinge- 
ment) and no external damping, 6, was applied to the blade, are re- 
ferred to as the "initial conditions. " A series of these initial conditions, 
one for each of three |i's chosen, was run for each of three chosen values 
of collective pitch, 0 q (see Table 2 for values). 

The procedure for obtaining the initial conditions at any ji for a 
given 0Q was as follows. The constants 0 for each location x were 
adjusted everywhere in the analog circuits to the appropriate value, 

9 = 6o + ®o^ (x-eo) 

Vq was adjusted in the circuits to the value determined from Equations 
(12) or (15). The thrust per blade, was measured for |j = 0. 0 

(V = 0.0), hover condition, by multiplying the time- averaged sum of the 
analog circuit signal outputs by the appropriate conversion constant. This 
measured thrust (see Table 2) was checked against the theoretical value 
determined from Equation (11), A close correspondence indicated all sys- 
tems were functioning normally. In addition to this check, the quantity 
0x^0 Vq Ax was measured from each of the five forcing circuits and 
summed to obtain an approximate torque per blade at hover, 

^1 2 

Q = Z 2 ® ^o) ^^i 

i= 1 

Note from Equation (13), that the drag terms were omitted in this measure- 
ment. This rather crude approximation vas utilized because of a hick of 
sufficient analog components. As a result, the measurement was made only 
for hover conditions. The measured value (see Table 2) was checked against 
the theoretical value at hover obtained from Equation (13). These checks 
assured that the behavior of the electromechanical model closely corre- 
sponded to the theoretical or math model, and that all systems of the 
electromechanical model were operating properly. 
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Initial conditions for tlie two forward flight configurations were ob- 
tained by adjusting the amplitude of V sin nt in the analog circuit to the 
desired value, then adjusting the d. c, term -Vttg + VQ to a level which 
brought Tg^^g back to the same value obtained at hover. While oper- 
ating the model, a response outimt near the tip of the blade, V9 (see 
Figure 5 for location), was analyzed simultaneously with an oscillo- 
scope and a spectrum analyzer. Components 0c* cyclic pitch 

amplitude, and i^q, the cyclic pitch phase, were adjusted in the circuit 
until the one-per-rev component of V9 was minimized, while at the same 
time trirnming the term Vq - VUg to maintain the desired thrust. It will 
be seen later that this procedure resulted in realistic values of the oper- 
ating conditions. When these initial conditions were obtained, scope 
trace photographs were taken of the velocity response V7 and V9 and the 
total a. c. component of the thrust I^ (see Figure 7). The latter com- 
ponent is the vector sum of the fluctuating part of the thrust at all blade 
increments. Note that in the expanded scale of Figure 7a, some one-per- 
rev component of V7 and V9 remains after minimization for V = 0, al- 
though no external impulses were being applied. This response results 
from small amounts of beam unbalance in the mechanical rotor. These 
records, along with the spectral plots of V9 to be described later, were 
used as a reference for the initial conditions and assured that we could 
return to those conditions at any time and Icnow the model was functioning 
consistently. 

With the model running under initial conditions, the values of 9j,, 
ijtQ, Tg^g, and -VUg + VQ could be measured. The value for Vq at the 
forward flight speeds was obtained from Equation (15) and the measured 
Tg^g. The value for Ug could then be determined from the measurement 
for Vq - V(Xg. All measured values for operation of the model in initial 
conditions are given in Table 2. Plots of thrust coefficient, Ciji/o, as 
a function of shaft tilt angle for prototype rotor blades operating at our 
values for and 9 q were obtained from W. H. Tanner’s work on rotary 
wing performance (Reference 11). Since the charts and graphs in Refer- 
ence 11 are referred to the axis of no feathering, it was necessary to 
account for the cyclic pitch in correlating these values to the quantities 
used in the present study, which are referred to the shaft axis. Essen- 
tially, this amounted to adding ti cyclic pitch to the rotor angle of attack 
(ttc, in the terminology of Reference 11) to obtain the shaft tilt angle 
shown in Figure 8. From the measured thrust, a "measured” thrust 
coefficient could be calculated for any initial condition using the equation: 



pc b R^Q^ 


Values for C«j</(7 obtained from measured thrust and the associated mea- 
sured cXg for any given initial condition are compared to the results of 
Reference 11 in the plot of Figure 8. The reasonable comparison of the 
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meaatired and theoretical values assured us that our model simulated a 
realistic prototype rotor blade operating under realistic conditions. 

Procedures for Test Runs 

Having established initial conditions as described in the previous 
section, the model was run for a variety of test conditions in -which a sig- 
nal simiilating vortex impingement was summed into the analog circuit. 
Furthermore, the effect of external damping on the blade response -under 
these conditions was investigated by feeding back the negative velocity re- 
sponse signal from certain points of the circuit to their corresponding 
exciter points. Thus, in effect, an equivalent viscous t-ype damper re- 
sulted. 


It was pointed out in a previous section that one would expect the 
electromechanical simulation to become invalid at some intermediate fre- 
quency because of the finite size of the forcing distribution increment. 

That is, for shorter bending wave lengths, both the simulated aerod-ynamic 
feedback and simulated damping feedback concepts break down. Some pre- 
liminary experimentation was performed to determine the useful upper fre- 
quency limit for the apparatus. It was foTUid that damping feedback had a 
strong tendency to amplify response above 10-per-rev, rather than provide 
positive damping. This result seemed plausible in view of the fact that the 
wavelength of the third bending mode is comparable to the force distribu- 
tion increment (see Figures 3 and 5). As a result, all responses fed back 
in the analog were low-pass filtered above 16-per-rev to reduce this ten- 
dency, It is estimated that only results below 10-per-rev are valid, how- 
ever. All results presented herein should be understood to include this 
restriction. 

Preliminary experimentation was also conducted to determine the 
general t-ype of responses that occurred for several different type of exci- 
tation pulses considered for demonstrating the effectiveness of the damping 
feedback. In addition to those described in the results to follow, a single 
pulse was also considered. However, for those p-ulse shapes utilized, the 
responses did not extend in time be-yond a single blade revolution, so that 
the results can be deduced from those of a pulse repeated at one-per-rev, 
and data on the single pulse was not acquired. 

Table 3 gives the matrix of test conditions investigated. Table 4 
explains in detail what the simulation function v(t) looked like in every 
instance. These different fimctions are coded in the matrix itself under 
a single "p” n-umber. Three levels of feedback were used. Values for 
the damping force per unit velocity, 6, were 0.0, 1.96 x 10^ 

N/m/sec, and 3.92 X 10^ N/m/sec, and are coded in Table 4 as 6 = 0, 

1, and 2, respectively. 
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Two different locatione for damping feedback were simulated. 

Outboard damping feedback was accomplished by amplifying the V9 velo- 
city signal to a desired level and feeding the result back into force 
exciter F5. Inboard damping feedback was accomplished by amplifying 
the V7 velocity signal to a desired level and feeding the restilt back into 
force exciter F4. 

The amplitudes and time widths of the pulses described in Table 4 
Were arbitrarily chosen to represent something like the inflow velocity 
disturbance caused by the impingement of a vortex shed from a preceding 
blade. The amplitudes were comparable to the steady inflow velocity Vq, 
and widths were of the order of 0. 6 meter or less in space. ^ 







SIMULATION RESULTS 


Preliminary Comments 

Results for a given collective pitch 0 q can conveniently be sepa- 
rated into three broad categories: steady state flight, the one-per-rev 
pulse condition, and the moving pulse condition. These categories relate 
to the type of function v(t) applied as defined in Table 4, and the results 
for each category will be presented separately below. Some escplanation 
of resTilts which tend to verify the validity of the simulation will also be 
presented. Further discussion related to results from the use of damping 
feedback will be presented in a following section, M all cases, the results 
are presented in terms of prototype parameters for greatest utility, al- 
though the measurements were, of course, performed on the model. 

Generally, the data is presented in several forms. Oscilloscope 
photographs of response time histories are given for key examples to show 
blade motion at several locations. Most of the data were acquired in the 
form of root mean square response amplitude and spectral decomposition 
of the response at V9. Recall again that V9 is the output of the second 
velocity transducer from the blade tip (see Figure 5). In some cases, 
oscilloscope time histories were also obtained for responses at V7, and 
are also presented. It would, of course, have been desirable to acquire 
data at all response points along the blade. However, this would have re- 
sulted in a further increase in an already large volume of data. Thus, it 
was considered feasible to form conclusions about blade response behavior 
from V9 results alone. 

All details of RMS and spectral response data are presented in the 
several tables of Appendix A. However, a substantial number of sample 
plots of these data for 0^ = 12° only, will follow below in order to provide 
a more convenient basis for discussion and forming conclusions. Root 
mean square amplitude data are presented in two forms: "as measured" 
(Table A-1), and "corrected" (Table A-2), narrow band results which are 
valid through 10-per-rev. These results were obtained by forming the com- 
bined RMS values from the first 10 components of the spectral data given 
in Appendix A, As described in previous sections, this range of frequency 
was considered valid for the analog simulation system. In addition, the 
corrected RMS response values include an accounting for certain crosstalk 
noise that appeared on V9 for those conditions where a pulse was applied 
(details to be. given later). All RMS amplitude response plots to follow are 
based on such narrow band results, and those for pulse applied conditions 
are based on corrected narrow band results. 
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Finally, spectral data were acquired up to l6-per-rev, and in 
some cases beyond this. Results are generally presented up to l6-per- 
rev, since this range was convenient for analysis purposes. However, 
those components above 10-per-rev should be considered invalid, and 
studied only for indication of how rapidly the simulation breaks down as 
frequencies increase. 


Steady State Flight 


Under steady state flight, no external pulse is applied to the blade, 
i, e., v(t) = 0. For this condition, the oscilloscope photographs of re- 
sponse time histories have already been shown in Figure 7. It is interesting 
to note that the response at both V7 and V9 is largely at Z-per-rev, and at 
|i = 0,42 the response at 2-per-rev is nearly equal in amplitude for both 
positions. By referring to Figures 3 and 5, this response can be concluded 
to occur principally from motion in the flapping and first bending modes, 
and the phasing must be such that the blade tip motion in each mode is out 
of phase, and the flapping response must be larger than than of the bending. 

Further data are shown in Figures 9 through 12. Figure 9a shows 
an example plot of the narrow band RMS amplitude response at V9 as a 
function of advance ratio |j, for fixed collective pitch Sg = 12®, and 
for outboard damping feedback set at several different levels. A com- 
plete set of these data for all flight conditions is presented in Table A-1. 
Further similar RMS response data are presented in Figure 9b for inboard 
damping. From these data it was decided that the effect of damping feed- 
back at the inboard location was negligible. As a result, hereafter damp- 
ing feedback was limited to the outboard location only. 


Figures 10 through 12 show example plots of the spectral response 
for V9 and steady state flight conditions. Similar data are tabulated in 
Table A-3 for all values of collective pitch. From these figures it can be 
seen that for the steady state conditions defined, the major components of 
response in the blade occur at the first four harmonics of blade rotation, 
i, e. , at 1, 2, 3, and 4-per-rev, From Figure 3 this can be seen to be 
response principally in the flapping and first bending mode. The only 
small to negligible responses at 6 and 7 -per rev, and near 11 and 1 2-per- 
rev, indicate that very little motion in the second and third bending mod|es 
respectively, is present in the response at V9, By inspection of Figure 5, 
it can be seen that V9 is located near a node for both the higher modes. 
Thus, motion in these modes is probably present in the overall blade re- 
sponse, but only a relatively small amount of it can be detected at V9. 

Even without the previously- described method of nulling the one- 
per-rev component for initial steady state conditions, one might again 
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■wonder why vibration at any frequency is present for the hover condition 
shown in Figure 10. As was previously explained, some vibration is ex- 
cited by residual unbalance in the mechanical rotor system. This vibra- 
tion is given by the results in Figure 10a. The amplitude of the component 
at one-per-rev would be far greater in Figures 11 and 12, which corre- 
spond to forward flight speed conditions, had the trimming not been used 
in the initial condition procedures. Note also that for forward flight 
speeds, the one-per-rav component could not be completely eliminated 
because the phase of the unbalance excitation does not coincide with that 
from forward flight excitation. 

Repeated Pulse 

For this series of tests, v(t) was implemented as described in 
Table 4 for p = 1 or 3, Essentially, these tests simulated a vortex of 
two different diameters and velocity amplitudes impinging at the F5 exciter 
location once every revolution of the blade. Figures 13a and 13b show 
time histories of the blade response at V7 and V9 for hover and for the 
higher forward speed, where no external damping has been included. The 
applied inflow velocity pulse in this case has an amplitude of 8.47 m/sec 
and a diameter of 0. 6 m. The affect of the inilow velocity pulse on the 
fluctuating component of the total blade thrust, I^;, can be obtained from 
the lower trace in each figure. Note, however, that inadvertently the 
pulse was lost off the bottom of the photograph in Figure 13b. For this 
case it wo'uld appear as a pulse downward with the same amplitude as in 
the bottom trace of Figure 13a. The upper trace in all these figures repre- 
sents only a timing reference signal for marking one-per-rev, at fit s 90**, 

Figures 13a and 13b should be compared directly with Figures 7a 
and 7c. Note again, however, that a different scale was used in Figure 7a, 

It can be seen that the effect of the pulse is to excite principally responses 
at several times the blade rotation, with smaller changes to the response 
in the range of one-per-rev. It is also apparent that a pronounced high fre- 
quency downward and upward spike occurs at V9 during the application of 
the pulse, followed by a ringing of the blade at frequencies of 3-per-rev and 
above. Only the ringing type response is apparent in the V7 signal. A 
careful scrutiny of the pulse part of the response at V9 led to the conclusion 
that this part of the signal is error due to cross-talk between the velocity 
response transducer at V9 and the two force exciters on either side of it 
which made up F5 (see Figure 5). Note that this essentially added noise did 
not extend all the way to V7. It was determined that the cross-talk could 
not be eliminated without significant changes in the apparatus. A numerical 
correction procedure was developed to eliminate approximately the effects 
of the cross-talk from the V9 RMS response data. First, a spectral analysis 
was performed for signals similar to V9 and It in Figure 13a, which 
approximated the presence of the pulse only. This spectral analysis had 
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components of rmiform amplitude at all integer -per- rev frequencies in the 
range of l6-per-rev and even higher. This uniform amplitude was con- 
sidered a noise level which was found to vary slightly with the type of pulse 
applied. Corrected narrow band RMS data were formed by subtracting out 
the effects of the noise from each component. The difference in \mcorrected 
and corrected RMS response values can be obtained by comparing corre- 
spqnding data from Tables A-1 and A- 2. The appropriate noise levels are 
noted on all spectral plots for conditions in which a pulse was applied. 

In retrospect, the above problem could also have been avoided by 
analyzing signals at V7 instead of those at V9. However, none of the final 
conclusions from the study would have been altered. 

Root mean square amplitude response for the p = 1 repeated pTilse 
condition and three levels of outboard damping appear in Figure 14, and 
corresponding spectral plots appear in Figures 15 through 17. The RMS 
data indicate that escternal damping decreases the response for the range 
through 10-per-rev. The spectral plots of Figures 15 and 17 also confirm 
the results presented in Figures 13a and 13b. That is, for a given level of 
external damping or suppression, the effect of the pulse is to excite re- 
sponses at frequencies of 3-per-rev and above, with smaller changes to 
the 1 to 2-per-rev components. 

Additional RMS amplitude and spectral data appear in Figures 18 
and 19 through 21, respectively for the p = 3 pulse condition. Basically, 
the same conclusions stated above apply also to these results with one ex- 
ception. It is apparent for this case that the higher feedback suppression 
level tends to increase the RMS level rather than decrease it at the lower 
speeds. From Figures 19 through 21, and Figure 3, it can be seen that 
this behavior is associated with a decrease in damping of the third bending 
mode near 12-per-rev. Apparently the phasing of the response of this 
mode is such that the damping adds energy rather than removes it. It 
should also be noted from Figure 5, however, that the F5 exciter, at which 
the external damping is applied, is located near a node for this mode. Thus, 
this result is in question, and probably is invalid because of the size of the 
finite blade increment as already described. In any event, this mode occurs 
above the 10-per-rev range, although its effects in this case are apparently 
being felt down into this frequency range. 

Moving Pulse 

For this series of tests, v(t) had the chax*acteristics described in 
Table 4 for p = 2 and 4. These tests were performed to measure the blade 
response for an impulse loading which intercepts a blade at a series of 
points traveling from the blade tip toward its root. The response time 
histories for V7 and V9 for hover and high forward speed with the p = 2 
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pulse condition applied are shown in Figure 22. In addition to this condition 
being formed by a repeated pulse » the polarity of the inflow impulse was 
also reversed. This can be seen from the upward protrusion of the pulses 
shown in the lower traces of Figures 22a and 22b, as compared with the 
corresponding downward protrusion of the pulse in Figure 13, for the 
p = 1 single repeated pulse. Note from the I|. trace that the inboard- 
most pulse was too weak to appear. 

It is apparent from these figures that a similar high frequency re- 
sponse is excited by this pulse condition when compared to the responses 
in Figures 13a and 13b, However, the responses at 1-2-per-rev also ap- 
pear to be influenced more, as might be expected. Further, the responses 
of the moving pialse appear to carry over into the next blade revolution, 
rather than damp out in a single revolution, as appeared to be the case in 
Figure 13 for the p = 1 pulse. 

The same comments made in the previous section about the cross 
coupling spike response also apply for this case. Therefore, a similar 
numerical correction was developed to eliminate approximately the effects 
of the cross-talk noise from the RMS response results. All plotted data 
are again based on corrected values. 

Root mean square amplitude of the response at V9 is shown in 
Figure 23 for the p = 2 condition. The results here are basically simi- 
lar to those shown in Figure 14 for the p = 1 condition. However, the 
levels at hover are higher in the moving pulse case. These results are 
confirmed by the spectral plots which appear in Figures 24 through 26. 

It is apparent that the increased vibration response at hover occurs in the 
1-3-per-rev region. Thus, these data support the similar conclusion 
formed from Figure 22, that the moving pulse has a greater influence on 
the lower frequency response than the single repeated pulse. 

A corresponding set of plotted data for RMS response and spectral 
decomposition are presented in Figures 27 through 30 for the p = 4 con- 
dition. Careful study of these data reveal a behavior similar to that for 
the p = 3 single repeated pulse condition shown in Figures 18 through 21, 
However, in going from p = 2 to p = 3, again the responses at lower 
speeds (including hover) appear to be influenced the most. This change in 
response is also concentrated principally in the 1-3-per-rev range, 

A response simulation for the p = 4 case could not be achieved 
for a collective pitch of 9 q = 11° (see Table A- 7 in the Appendix). For 
this case, the system went unstable at the high forward speed ([i = 0,42). 
As a result, data for this case were taken for the p = 5 condition defined 
in Table 3. The system remained stable for the latter case. It should be 
emphasized that this instability moat likely resulted from the improper 
simulation above 10-per-rev, rather than any response in the valid fre- 
quency range. 
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DISCUSSION 


The preceding sections have already included a substantial amount 
of discussion of results presented. However, that discussion has been 
aimed principally at pointing out aspects of the data which tend to verify 
the validity of the simulation, and to show that various forms of results 
obtained were consistent. The purpose of this section is to emphasize the 
utility of the results presented in this report for consideration in the 
potential design of active feedback suppression systems in actual rotary 
wing aircraft. Of necessity, some of the discussion will expand upon that 
already presented. 

Generally, the results show that the use of an active response 
suppression device can be quite effective in reducing the overall response 
of the vibrating rotor blade. However, the location of the device obviously 
has a strong influence on those modes that can be suppressed. It is appar- 
ent that the device must be located near the tip of the blade to be the most 
effective. This result could have been deduced intuitively; however, one 
might not h ^ve realized how dramatic the effect was in moving from 
X = 4. 42 meters out to x = 5, 6 meters, as shown in Figure 9. It woxild 
even be more desirable to select a position even nearer the blade tip, 
which woTild correspond to an antinode position for most of the lower modes. 
Note that the F5 position in the present study nearly coincides with nodes 
for both the second and third bending modes. As a result, little suppression 
of these modes could be accomplished. 

After a careful study of all of the data near the end of the program, 
it became apparent that the early conclusion about the overall ineffective- 
ness of inboard damping probably was premature. The outboard damping 
was effective in suppressing the 1-3-per-rev flight loads response which 
was primarily blade flapping and first bending. The outboard feedback was 
active in both of these modes. However, the situation appeared to be re- 
versed for the impulse suppression. The outboard feedback was not active 
in the second and third bending modes as a result of the feedback location 
being at the node points of these modes. Therefore, the outboard feedback 
location tended to be ineffective in suppressing the higher frequency blade 
response associated with impulse loading. Unfortunately, the results are 
inconclusive regarding inboard feedback location suppression of impxilse 
response. However, the inboard location had higher participation in the 
second and third modes and, therefore, could be expected to be more effec- 
tive in impulse response suppression. 

In most of the spectral data presented, it appears that the 1 and 3- 
per-rev components were significantly reduced by the outboard damping. 
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but the 2-per-rev component was essentially unaltered for most cases. 

The reason for this is not clear. It has been mentioned that the 2-per-rev 
component must be comprised principally of the flapping and first bending 
modes. However, these occur very nearly at 1-per-rev and 3-per-rev, 
respectively. Thus, if these components are reduced by damping, it seems 
reasonable that the 2-per-rev component should correspondingly reduce. H 
is plausible that the addition of e^cternal damping causes a shift in the phas- 
ing of the responses of the two modes such that the vector stim of the two 
still is affected relatively little, although the individual amplitudes are re- 
duced. Interaction with the effects of the cyclic pitch may also be the 
source of this behavior. In the present case, cyclic pitch was adjusted for 
a minimum 1-per-rev initial condition and then e3d:ernal damping was applied. 
Further adjustment of cyclic pitch thereafter may have reduced tiie 2-per-rev 
component, had this been attempted. In any event, this result can be most 
significant in design. 

Another lesson must be learned from the effects of the location of 
the suppression device on damping characteristics. Considerable care 
had to be exercised with the model rotor blade to filter out higher mode 
response above l6-per-rev to avoid instability in those modes. This 
difficulty was caused by the use of a finite increment in representing the 
aerodynamic and damping force on the blade. Thus, in practice, care 
will have to be taken to assure that the response sensor and reaction force 
applicator are located very near oach other in order to avoid a similar 
instability from occurring in a prototype system. 

In this present study, the diameter of a vortex was assumed to be 
about 0.25 meter or smaller. It is apparent from the results that such a 
narrow repeated rapping of the rotor caused principally high frequency 
response, although the moving pulse did have a relatively significant ef- 
fect on responses below 4-per-rev. A careful investigation should be 
pursued to determine whether larger diameter vortices are realistic. If 
they are, it is apparent that considerably greater lower frequency re- 
sponses will be excited by such larger diameter (wider pulse vortices). 

At the same time, an investigation of realistic vortex amplitude is also 
warranted, since this parameter has a significant influence on the re- 
sponse results. 

The results of the present study can be used to provide some idea 
of design parameters for external active modal suppression. The RMS 
results can be used directly to arrive at a damping effectiveness coeffi- 
cient in terms of percent vibration reduction per unit of damping force 
(N/m/sec) applied. From the several RMS plots it is obvious that this 
coefficient would vary with |a as well as the type of pulse excitation con- 
sidered. However, an average value would be of the order of 30/1,96 X 10^ 
percent per K/m/sec, or about 0. 015% per N/m/sec applied. 


Z6 


Various limitations of the present apparatus have been pointed out 
in the discussion throughout this report. However, it has still provided 
useful information directly applicable to design of actual active modal sup- 
pression devices for helicopter rotor blades. Having the benefit of this 
first experience in studying this problem, a variety of further studies sug- 
gest themselves. Obviously the significance of higher frequency response 
eind possible instability can be investigated by using more of the total of 
nine force increments available with the mechanical blade. Unfortunately, 
this dictates the use of a more elaborate analog computer system. The 
possibility of using more sophisticated aerodynamic theory can be con- 
sidered. Considerably more data could be taken and analyzed from other 
blade response positions, even with the present facility. It is particularly 
obvious now that more work with inboard damping and impulsive loading 
conditions would yield further useful results. The matter of avoiding 
cross-talk between force and response channels could be better effected. 
With a little more effort, a vortex impingement pattern appropriate for the 
prototype rotor at |i = 0. Z6 could be achieved. A more elaborate analog 
setup would also have allowed a more accurate confirmation of total torque 
as well as total thrust. Even so, it is felt that the effort has allowed a 
preliminary, but significant, demonstration of the effectiveness of active 
feedback suppression by means of a rather unique technique. 


TABLE 1(a). MODEL PARAMETERS 


n /n =2.5 

m p 

V /V = 0.216 

m p 

Bending Ntimber 

El 

Flapping Inertia 

3p pcjj 
m 


Pm/Pp = 1.428 

Fj^/Fp = 4.98X10^'^ 

0.00923 

2. 73 


TABLE 1(b). helicopter FIXED PARAMETERS 


Blade Length, S, 

Ntumber of Blades, b 
Flap Hinge Eccentricity, 6 q 
B lade Chord, c 
Air Density, p 

Mas s per Unit Length, m 
Blade Stiffness, El 

Rotational Speed, O 

Max, Blade Twist, 0j 

Slope of Drag Coefficient, Cjjq 
(B lade NACA 0012) 

Blade Tip Speed, (eo+A)fi 

Slope of Lift Coefficient, C^q 
(B lade NACA 0012) 


Prototype 
6.01m (19.7 ft) 

4 

0.354m(l.l6ft) 

0.533m(1.75ft) 

1.23 kg/m^ 

(1.15 X lO-’^lb-sec^/in^) 

13.6 kg/m (9 lb /ft) 

1.58 XlO^N-m^ 
(5.5xl0'^lb£-in2) 

300 RPM 

- 8 ° 

.092/rad. 

200 m/sec (655 ft/ sec) 
6.13/rad. 


Model 

0.619rti(1.70ft) 

4 

3.06 X lO-^m (0.100 ft) 

4.61 X 10“2m (0.151 ft) 

1.75 kg/m^ 

(1.640 Xl0"'^lb-sec2/in'^) 

0.145 kg/ra (9.60 Xl0“2lb/ft) 

0.588 N-m^ 

(2.05 lb£-in2) 

750 RPM 
0.092/rad 

43,2 m/sec (141 ft/sec) 
6.13/rad 
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TABLE 2. MEASURED PROTOTYPE HELICOPTER VARIABLE PARAMETERS 


Collective Pitch, 9q, deg 









( at blade root) 

10 



11 



12 


Advance Ratio, u 

0 

0.26 

0.42 

0 

0.26 

0.42 

0 

0.26 

Forward Speed, V, m/sec 

0 

52.4 

84.1 

0 

52.4 

84.1 

0 

52.4 

(ft/ sec) 

(0) 

(172) 

(276) 

(0) 

(172) 

(276) 

(0) 

(172) 

Shaft Tilt Angle, ttg, deg 

0 

-6.10 

-3.72 

0 

-6.75 

-3.76 

0 

-6.36 

Cyclic Pitch Amplitude, 
0c, deg 

0 

2.80 

2,56 

0 

3.57 

3.33 

0 

4.43 

Cyclic Pitch Phase, 
<3eg 

1 

0 

3.95 

4.65 

0 

3.25 

3,25 

0 

4.00 

Inflow Velocity, Vq, m/ sec 

6,75 

0.875 

0.545 

7.97 

1.23 

1.08 

9.13 

2.31 

(ft/ sec) 

(22.15) 

(2.87) 

(1.79) 

(26.16) 

(4.04) 

(3.56) 

(29,95) 

(7.57) 

Thrust Per Blade, Tg, 









N X 10-3 

4.09 

4.09 

4.09 

5.19 

5.19 

5.19 

6.62 

6,62 

(Ibf) 

(920) 

(920) 

(920) 

(1166) 

(1166) 

(1166) 

(1488) 

(1488) 

Torque Per Blade, Q, 









N-m X 10-3 

2.44 



3.62 



4.92 


(ft-lb£ xlO-3) 

(1.80) 


-- 

(2.67) 

-- 

-- 

(3.63) 

-- 




TABLE 3. MATRIX OF TEST CONDITIONS 


4° U p* 6* 

deg 

f, 

deg P ® 

®0 c 

dsg P ‘ 

12 0 0 0, 1, 2 

11 0 0 0, 1, 2 

10 0 0 0, 1, 2 

1 ” 

1 

1 ’• 

2 '» 

2 

2 " 

3 " 

3 " 

3 " 

4 ” 

4 " 

4 " 

0.26 0 ” 

0.26 0 ” 

0.26 0 " 

1 " 

1 " 

1 " 

2 " 

2 " 

2 *' 

3 " 

3 " 

3 " 

4 " 

4 " 

4 " 

0.42 0 " 

0.42 0 " 

0.42 0 " 

1 " 

1 ” 

1 " 

2 " 

2 ” 

2 " 

3 

3 ” 

3 •' 

4 ” 

4 " 

5 " 

* See Key to Test Matrix on next page. 





TABLE 4. CHARACTERISTICS OF APPLIED 
PULSES AND GUST SUPPRESSION 


Applied Pulse V(t) 


Pulse Pulse 

Condition Amplitude 


Pulse 

Width 

millisec 


Point and Rate of Application 


8.47 

(27.8) 

8.47 
(27. 8} 

25.4 

(83.4) 

25.4 

(83.4) 

16.9 

(55.6) 


No pulse applied 

R-6o = 0.88 (at F5), and 1-per-rev* 

R-eo = 0.88, 0.69, 0.50, 0.36, 0,21'^ 
(at F5, F4, F3, F2, FI) & 1 each per rev 

R-Oq = 0.88 (at F5), and 1-per-rev 

R-eo = 0.88, 0.69, 0.50, 0.36, 0.21 
(at F5, F4, F3, F2, FI) & 1 each per rev 

R-eo = 0.88, 0.69, 0,50, 0.36, 0.21 
(at F5, F4, F3, F2, FI) & 1 each per rev 


Gust Suppression Feedback 


Damping Damping 


Condition 


Constant 

N/m/sec 


1.96 X 103 
3.92 X 103 


"Since the response of the blade is unaffected by the applied pulse for times 
longer than one blade revolution (see Figure 13 for example), a single 
applied pulse produces the same response in the blade as appears for one 
revolution of the p = 1 condition, 
t See Figure 4 for source of these data. 




















{ a ) Rotation component in place of blade 



( b } Span-wise velocity component perpendicular to blade 


Figure 2. Relationships between angular veixity or span -wise velocity 

and the blade deflections 


lO 

ROTATIONAL FREQUENCY, ( rpm ) 36 . 

Figure 3, Natural Frequency Diagram of Model Helicopter 

Rotor Blade 






POINT ON BLADE INTERCEPTED BY VORTEX 



Figure 4. Intercept Point of Vortex on a Following Blade 





1st MODE 


BLADE. X 


2nd MODE 

3 rd MODE 
V5 V7 

I V6 I V8 



L. 

L F2 F3 

F- 


-ru uu I u 80 90 

PERCENTAGE OF DISTANCE FROM ROTOR AXIS. (%) 


EXCITER 
POINTS 


Blade IVlode Shapes and Location of Exciter and Response Points on Blade 





^ V— 


-J V. 


SOURCE OF ANALOG ANALOG INPUT 
INPUT SIGNALS SIGNALS 


ANALOG CIRCUIT ELEMENTS 


V- 


OUTPUT TO F5 
ON BLADE 


OUTPUT TO BLADE 
EXCITER 


Figure 6. Block Diagram of Apparatus to Generate Excitation at F5 
















, Timing Pulse ^at Ot = 90 
^ ( 1 per Rev ) 


[ 0.18 m/sec,V7 
0.18 m/sec,V9 

J 1.02 X 10' Newtons, 1^ 



I 2.56 X 10^ Newtons, 1^ 


J 1.81 m/sec,V7 
, 0.90 m / sec, V9 



Timing Pulse 


’[ 1.81 m/sec, V7 
' 1.81 m/sec, V9 
J 5.12 X 10^ Newtons, 1^ 


NOTE: V7 and V9 have different scale in Figure 7a 
and Figures 7b and 7c 


igure 7. Response of Blade for Initial Conditions of Hover and 
Forward Flight for = 300 r pm u(t) =O0o=12“, 
No External Damping 





THRUST COEFFICIENT, Cjlo 



Figure 8. Comparison of 


MEASURED VALUES 
o 00 = 12“ p=0.26 

• =12“ =0.42 - 

□ =*lf =0.26 

■ =lf =0.42 

O =10" =0.26 _ 

♦ =10" =0.42 



! I I I 

6 7 8 9 1 

Erctg, degrees 


and Modeled Thrust Coefficient 



VELOCITY, V9. ( m/sec } RMS VELOCITY, V9. (m/sBC ) RMS 



( a ) Outboard Damping Feedback, 6, x » 5,6 meters 



(b) Inboard Damping Feedback, 6, x- 4.42 meters 
Figure 9. Blade V9 Response for 0 g ■ 12*. v(t) » G, p - 0 
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ROTATIONAL FREQUENCY, per Rev, ( 1 per Rev - 5.0 Hz ) 


Figure 10. Spectral Plots of V9 Response, M *0, u(t) -0, 0o“ 12* p = 0 











Timing Pulse at Qt - 90° 
( 1 per Rev ) p 



1.81 m/sec,V7 


1.81 m/sec, V9 


5.12 X 10^ Newtons, 
p *= 1, Pulse Applied 


It 


Timing Puise 
I ( 1 per Rev ) \ 


(b) m-0.42 



1.81 m/sec,V7 
1.81 m/sec,V9 

5.12 X 10^ Newtons, 


Figure 13. Response of Blade V7, V9 for Q » 300 rpm 

v(t) = 1 per Rev Pulse, = 12°, No External Damping 
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ADVANCE RATIO, u 


Figure 14. Blade V9 Response for Outboard Damping Gq » 12**, v(t) * 1 per Rev Pulse ( p * 1 ) 







ROTATIONAL FREQUENCY, per Rev, ( 1 per Rev - 5.0 Hz ) 


-40 NOISE 
LEVEL 


Figure 15. Spectral Plot of V9 Response, u = 0, v (t) ■ I per Rev, 0o = 12” p’l 









ADVANCE 


Figure 18, Blade V9 Resposne for Outboard Dampl 
















( a ) 


Timing Pulse at Qt - 90' 
( 1 per Rev ) ; 



^ 1.81 m/sec,V7 
^ 1.81 m/sec, V9 


Timing Pulse 


5.12 X 10^ Newtons, Ij 
p - 2, Pulses Applied 


( 1 per Rev ) r 


(b) M '0.42 



1.81 m/sec,V7 
J 1.81 m/sec, V9 

21 5.12 X 10^ Newtons, 4 


Figure 22. Response of Blade : V7, V9 versus Time for Q =300 rpm 
\) (t) “Moving Pulse, eQ = 12®, No External Damping 



VELOCITY, V9, (m/s«c)RMS 












Ln 

OD 



Figure 27. Blade V9 Response for Outboard 



ing 0 Q = 12®, v(t) » Moving Pulse ( p = 4 ) 
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-30 

-50 

-60 
“20 

-30 
-40 
-50 
-60 
-20 
-30 
-40 
-50 
-60 

0 

ROTATIOWAL REQUENCY, per Rev, ( 1 per Rev - 5.0 Hz ) 



Figure 30. Spectral Plots of V9 Response, p * 0.42, u(l | - moving pulse. 0 , 
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TABLE A- 1. narrow BAND RMS AMPLITUDE RESPONSE AT V9 
(Frequency Range 1 - lO-per-rev) 





















TABLE A-2. NARROW BAND CORRECTED RMS 
AMPLITUDE RESPONSE AT V9 

(Frequency Range 1 - lO-per-rev) 


0o = 12° 


9o = 11" 


8o = 


0.26 


0.42 


N/m/ 

i sec X 10"^ 


0 . 0 
1. 96 
3.92 
0 . 0 
1.96 
3.92 
0.0 
1.96 
3.92 


0. 0 I 


0.42 


0.26 


0.42 


0.0 
1.96 
3. 92 
0 . 0 
1.96 
3. 92 
0.0 
1.96 
3. 92 


0 . 0 
1.96 
3.92 
0,0 
1.96 
3.92 
0.0 
1.96 
3. 92 


m/sec m/sec m/sec m/sec 


P_=_P. 
0. 068 
0.019 
0. 009 
0. 317 
0. 192 
0. 132 
0. 782 
0.490 
0. 315 


E^_l 

0. 065 
0. 025 
0 . 021 
0. 286 
0. 182 
0. 117 


0. 287 


0. 029 
0. 381 


0. 113 


0. 514 


0.213 
0 . 060 
0. 185 
0. 367 
0 . 200 
0.227 
0.759 
0. 562 
0.404 

;0. 256 

0.060 

jo. 175 

! 0.422 
0. 180 
0. 135 
0 . 862 
0.483 
0. 361 


0. 336 


.228 


- 201 
. 047 
. 177 
.271 
. 135 
, 112 
. 576 
. 324 
.291 


0.220 
0. 052 
0. 105 
0. 302 
0. 166 
0. 135 
0.495 
0 . 298 
0. 231 


0. 240 
0 . 061 
0. 160 
0. 358 
0. 148 
0. 173 
0.632 
0. 332 
0. 312 


0.285 
0.058 
0. 135 
0. 933 
0. 300 
0, 169 
0. 572 
0.277 
0. 187 



P = .0. 

0. 

094 

0. 

031 

0. 

017 

0. 

175 

0. 

134 


0.262 
0. 096 
0. 122 
0.541 
0.212 
0. 171 

p = 2 


V9 

m/sec 


0. 226 
0. 076 
0. 159 
0.476 
0. 115 
0 . 162 
0.573 
0.233 
0.217 























TABLE A- 3 


V9 SPECTRAL RESPONSE AMPLITUDES 
For v<t) = 0, p = 0 



bb 
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table A-4. V9 spectral response amplitudes 

For v(t) = 1 Por Rev, Amplitude 8.47 m/Sec, p = I 


u = 

mm 

0.0 

o 

o 

0.26 

tm 

IQQH 

0.42 

0*42 

0,42 

6(N/m/ 
aec it 10'^ = 

0 

1.96 

3, 92 

H 

1.96 

3.92 

■ 

1.96 

3*92 

e„ . 12“ 

(0 dB =: 

1. 81 m/aec rma) 







Rot.Freq* 
(per Rev*) 

-dB 

mm 

-dB 

-dB 

-dB 

-dB 

-dB 

-dB 

-dB 

1 

28.0 

42.0 

46.0 

24.6 

35.2 

40.8 

14.2 

26,5 

32*2 

Z 

36*5 

38.7 

43.2 

16.0 

19.2 

22.6 

8.5 

11*4 

15.5 

3 

31*5 

40,0 

43.5 

33*0 

62.3 

48. 6 

21.7 

34,8 

44.4 

4 

33,8 

36,2 

40. 0 

32.0 

35.8 

40*5 

33.6 

40.1 

44,0 

5 

39.0 

39.2 

42.2 

37.2 

39.0 

43.0 

36.8 

37.8 

42*8 

& 

40.0 

37.8 

39. 5 

37.2 

37*2 

40.8 

36.2 

37.8 

40*7 

7 

38.0 

37.0 

39.8 

37.2 

38.0 

41.3 

38,3 

41*0 

47.5 

8 

39v0 

37*0 

37.2 

37.2 

36. 3 

38.9 

37.2 

39.0 

40*8 

9 

38*8 

35,0 

35.0 

3 7*2 

36.0 

37, 3 

36.8 

38.3 

39.6 

10 

40.0 

35. 1 

35.2 

37.2 

35,0 

35.9 

38.2 

37.3 

37,7 

11 

39.2 

35.3 

34.0 

36.5 

35*0 

34. 1 

37,2 

37.2 

36.0 

12 

39,5 

35.8 

35, 1 

37. 0 

34.8 

35. 0 

35,4 

36.0 

37.7 

13 

39.5 

34,2 

30.2 

37.0 

33, 5 

31.3 

36.3 

37.0 

34.4 

14 

38.5 

33.8 

26.5 

36*2 

34.0 

28.8 

37.8 

36.2 

31.5 

15 

38.5 

35.0 

26. 0 

38* 1 

34.3 

27. 8 

38. 8 

37. 5 

31.5 

16 

39.0 

34.3 

25. 1 

37. 3 

34. 8 

26.2 

37.2 

38.2 

30.6 

9o = 11“ 

1 

18.5 

31.3 

34.0 

20. 6 

33.4 


23.5 

39.5 

40. 3 

2 

33*2 

36. 5 

37.7 

20*8 

22.5 


13.0 

15.0 

18.6 

3 

31.0 

39.0 

41. 0 

19.7 

30*2 

35. 1 

12.5 

23. 8 

19.5 

4 

36.2 

38.0 

41.0 

36. L 

35.4 

40.2 

29.0 

33.3 

36, 3 

5 

36.8 

35.0 

35.3 

34.5 

33.5 

37,2 

33. 0 

34*5 

38.0 

6 

39.2 

36.5 

33.6 

36.9 

34.7 

35.7 

37,0 

34.8 

36.0 

7 

38*0 

34.2 

31. 1 

37.3 

33.8 

34.8 

36.0 

32,5 

34*5 

6 

40.5 

34.9 

29.4 

37.6 

33.3 

32,0 

37.2 

33.0 

31.5 

9 

39.2 

35.0 

26.5 

38.7 

33.0 

37.7 

37.0 

32.8 

27, 1 

10 

40.2 

35*0 

20.8 

37.6 

33.0 

23. 3 

37.5 

32,5 

23,2 

11 

40.6 

35*5 

21.9 

39. 1 

33.7 

23,8 

38.2 

33,0 

23,5 

12 

41,0 

36.0 

23.0 

36. 3 

33.5 

22,6 

36.0 

32.0 

18.6 

13 

40.1 

35.5 

22.3 

38.7 

33.2 

24.5 

37,5 

32.5 

24.2 

14 

40.5 

37.0 

27.0 

39.5 

35,0 

29.0 

38.2 

35.0 

27.0 

IB 

40.8 

37.6 

30*4 

39.2 

34.7 

31.0 

38,8 

34.0 

29,8 

l6 

41.5 

38.3 

33.0 

39*5 

37.7 

33.0 

38.2 

34.5 

31.5 

0. = 10“ 






1 

20.3 

33.7 

39.0 

18.9 

32.2 

36,6 

15,7 

30,8 

33,3 

2 

33,0 

' 37,0 

40.2 

24,5 

' 26.0 

! 28.9 

17, 2 

! 19.2 

22,7 

3 

31*5 

1 39.8 

42.5 

20,2 

30.5 

35.5 

12*8 

24*3 

29* 3 

4 

35*9 

1 36.8 

41.6 

36. 0 

37.7 

39. 2 

28,0 

32.5 

35.2 

5 

37*8 

35.7 

39.3 

36* 1 

34.7 

37.0 

33*6 

35.5 

37* 1 

6 

39,1 

37.0 

38.3 

37,6 

34. 7 

35.6 

36*7 

34,5 

35,0 

7 

38.4 

34,3 

36. 0 

36.8 

33.2 

33. 8 

36.1 

32. 7 

33.0 

8 

41.0 

35.7 

34. 1 

38.8 

35.0 

32.6 

37.5 

34.0 

31.5 

9 

40.6 

34.1 

29.0 

39.0 

33,5 

27- 5 

38* 1 

33.2 

26.7 

10 

41,0 

35,6 

24.8 : 

39,2 

34,0 

24.6 

39-0 

33,7 

24.7 

11 

40.3 

36.3 

25.0 

40.5 

35,2 

26.0 

38.0 

34.3 

27.2 

12 

41.0 

36*5 

26*8 ' 

38*0 

35.6 

27.4 

35,5 

32*5 

27*5 

13 

40.8 

36.6 

27*0 

39*3 

34.4 

27.5 

38.3 

34.4 

30*0 

14 

41. 1 

37.3 

32,2 ! 

37*0 

35,4 

31. 3 

39.3 

34.5 

34*4 

15 

41.8 

37*4 

33,8 

40*0 

36.0 ! 

34. 1 

38,1 

35.3 

37.0 

16 

42.5 

1 

38*5 

36.3 

38. 9 

1 

37.3 

1 

36.5 

38*6 

36.5 

i 

39.2 


















































TABLE A-5. V9 SPECTRAL RESPONSE AMPLH’UDES 
For v(t) - 1 Per Rev Pulae, Amplitude 25,4 m/Sec, 

3 


Bo = ’2 


Rot.Freq. 
(per Rev.) 


(0 dB = 1.81 m/eec rtns) 


8 

35.0 

9 

35.7 

iO 

35.5 

LI 

36.0 

12 

35.8 

13 

36.0 

14 

36.0 

15 

36.9 

L6 

36.8 

O- 11“ 


-dB 

29. 

0 

14. 

0 

26. 

2 

28. 

4 

29. 

6 

27. 

5 

26. 

0 

25. 

3 

22. 

7 

21. 

0 

25. 

7 

21. 

0 

35. 

0 


!5. 

0 

14. 

5 

18, 

3 

31, 

4 

31. 

4 

33. 

0 

33. 

5 

34. 

4 

34. 

1 

34. 

5 

35. 

2 

33. 

6 

34. 

7 

34. 

4 

35. 

7 

35. 

8 


29. 

1 

25. 

7 

29. 

0 

32. 

3 

30. 

7 

31. 

0 

o 

cn 

2 

30. 

0 

30. 

0 

29. 

8 

30. 

7 

30. 

I 

30. 

1 

31. 

.5 

32. 

, 1 

33. 

0 


31. 

8 

28. 

9 

31. 

3 

33. 

2 

31. 

5 

30. 

7 

29. 

5 

26. 

7 

34. 

0 

21. 

2 

23. 

6 

25. 

0 

29. 

4 

36. 

0 

38. 

.0 

37. 

.3 


14.3 

28.0 

IB.O 

18.5 

12. 1 

23. 1 

27.4 

30.2 

30.7 

31.2 

34,0 

31.8 

34.4 

30.0 

34.4 

29.7 

34.9 

30, 1 

35.0 

30.0 

35,3 

31.0 

33. 1 

30. 1 

34. 1 

30.0 

35.5 

31.6 

35.0 

32.0 

36.0 

33, 1 


30. 

1 

22, 

0 

27. 

0 

32. 

0 

30. 

5 

30. 

3 

29, 

3 

26. 

3 

24, 

5 

21. 

7 

26, 

2 

26. 

0 

33, 

6 

39. 

2 

35. 

9 

35. 

,0 



















































TABLE A-6. V9 SPECTRAL RESPONSE AMPLITUDES 
Tor v{t) s Moving Pulse, Amplitude 8*47 m/Sec, 


6(N/m/ 
sec X 10“^ = 


So 12° 


Rot*Freq. 
(per Rev*) 


(0 dB = 1*81 m/sec rms) 


0.26 

0.42 

0,42 

0*42 

3.92 

0 

1.96 

3.92 


-dB 

-dB 

17* 5 

30. 6 

17*6 

20.6 

17*8 

29.4 

34,7 

36.0 

36*0 

35. 3 

37* 1 

37.5 

36, 9 

35.2 

38, 0 

35, 6 

37.8 

34, 3 

37*7 

34. 3 

1 38* 1 

33. 3 

' 38. 8 

34.0 

1 37* 8 : 

34. I 

37. 0 1 

32.4 

; 36*5 

32.0 

37.2 

33.0 


-dB 

-dB 

27* 3 

29* 1 

12*6 

16*0 

15.0 

27.2 

31.5 

31*0 

35. 0 

33, 1 

33.6 

32* 3 

31.5 

30.0 

1 34.0 

31*0 

33*7 

29. 1 

33.8 

27*5 

35,0 

29*4 

30* 7 

27. 2 

32, 6 ! 

28*7 

32. 5 

27. 3 

32.8 

25.6 

34*0 

29*4 


18, 0 

31. 0 

35. 8 

19*2 

32*3 

37. 1 

22*5 

36.4 

31*5 

34.4 

38, 0 

22*5 

24* 0 

27* 3 

23. 6 

16.0 

27. 5 

35*8 

39. 8 

18*0 

28. 0 

33. 3 

11*5 

22.6 

37. 0 

37.3 

42.0 

34* 1 

26, 5 

40. 1 

26.0 

32*4 

37.0 

35.4 

37.9 

34*2 

24.3 

37, 3 

33.0 

36.0 

41,0 

37. 0 

39.4 

39. Z 

36. 0 

37*0 

37*5 

35.5 

37. 3 

34.0 

35,2 

36*8 

33*3 

35*0 

36.8 

33.5 

39. 7 

35*0 

33. 2 

38.4 

34. 6 

32,6 

38.0 

32*9 

40. 3 

34.9 

29. 7 

38. 0 

33* 3 

28.8 

37.6 

32*7 

41,0 

34. 3 

24. 3 

38,0 

33.6 

23.4 

37*5 

33.0 

39.8 

34.8 

23, 5 

38*7 

33*5 

22* 5 

37* 8 

32,9 

38,0 

33.4 

24.0 

35.2 

30*5 

20. 1 

34*0 

29.0 

40. 2 

35*5 

26. 0 

37. 1 

34. 5 

25.0 

29*2 

33*7 

40.6 

37.0 

30. 5 

37*0 

34.3 ' 

27.9 

37,0 

33*8 

40. 1 

36. 0 ' 

33.0 

36* 1 

34. 5 ■ 

30*3 ' 

37*6 

33.7 

40,2 

37.5 

35. 0 

37.0 ' 

36.0 

32. ? 

38* 0 

35*0 


33.5 

37.2 

16* 7 

30*5 

3^8 

14*7 

29.2 

32*0 

35*1 

38.8 

26.0 

27.6 

3U3 

IB* 1 

19.9 

23*2 

36,0 

40. 6 

18. 8 

28.4 

33*0 

12.0 

22*5 

27*0 

37, 5 

41,0 

33.0 

35, 0 

38. 0 

26.0 

31.7 

35.0 

36.2 

39.6 

35. 1 

35* 1 

37.5 

33*7 

35.6 

35*5 

39.3 

39. 1 

38*2 

3b. 5 

36. 7 

38,4 

35.2 

34.8 

34.6 

35*7 

36* 1 

33. *7 

33.5 

36*4 

32*7 

32. 1 

36,0 

33.4 

38.6 

34*6 

32.2 

38.5 

33*7 

31.0 

35, 0 

29.2 

39.0 

33*5 

27.5 

37*9 

32.3 

26*4 

35.4 

25.6 

38.8 

34. 1 

25*7 1 

39.0 

33.0 

25*5 

36. 0 

25*0 

39.2 

35. 3 

25. 8 

39.5 

34.0 

26.8 

34* 0 

25. 0 ! 

33. 2 

29. 7 

22, 2 

32*2 

28.7 

21, 3 

38.2 

27.8 1 

38, 3 

35.4 1 

30.0 

40. 0 

34.5 

32.6 

38*0 

32.7 

39.5 

36.5 ! 

32*2 ' 

38*0 1 

34.6 

37* 1 

37.5 

34.7 , 

40. n 

35*6 

34*5 

39.0 

35*5 

39*4 

38*3 

35*5 

30.0 

36, 6 , 

36* 8 

37*6 

35.0 

39* 1 

































TABLE A-7. V9 SPECTXlAL RESPONSE AMPLITUDES 
For v(t) Moving Pulso, Amplitude 2B*4 m/Sec, 

. 5 



0O = 


(0 dB = 1,81 m/scc rms) 


0. 26 

0,42 

0.42 

0.42 

3.92 

0 

1.96 

3.92 


-dB 

-dB 

-dB 

-dB 

38. 6 

15*8 

25.7 

27.0 

24.9 

9.5 

11.6 

14.6 

30,4 

9.5 

21* 1 

24.8 

34.2 

25,0 

27,0 

28.8 

34,8 

30.3 

30.2 

30. 1 

33*9 

32.8 

29.0 

28.6 

29,8 

32.8 

27.8 

26. 1 

29.4 

33.6 

28.0 

25*8 

24.0 

33.8 

27*8 

23*4 

27. 0 

34.8 

29,9 

24.5 

31.5 

35.5 , 

30.6 

30,6 

26.4 

30.7 * 

24,3 

18.1 

32. 3 

35.0 

29.7 1 

33* 1 

34.9 

34.3 

30.0 ! 

35.0 

34.6 

33.5 ' 

30.8 

34.0 

30.6 

33,2 

30.5 

33*0 
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SYMBOLS 


b 




E 



I 

Jt, 

m 

P 

^total 

Qx 

R 

t 

T 

T 

J-avg 

U 

Up 

Ut 


number of blades 
blade chord , m 

slope of drag -coefficient versus a and lift- coefficient 
versus a curves, respectively 

spanwise distance from rotation axis to blade root 
(eccentricity), m 

blade elastic modulus, N/m^ 

eo/E- 

airfoil section drag and lift, N/m. 

blade spanwise bending moment of inertia, 

blade span, m 

blade mass per tuiit length, kg/m 

identifying number for applied v(t) function 

total aerodynamic torque about rotation axis, N-m 

airfoil section aerodynamic torque about rotation axis, 
N-m/m 

A + 6q , m 

time, sec 

airfoil section thrust, Fjj cosaj. - Fjj sinttj., N/m 
Total thrust averaged over one revolution, N 
relative wind, ^U>p^ + Up^ , m/sec 

total upward velocity in Yq direction (Figure 1), m/sec 

total velocity, positive in negative Zq direction (Figure 1), 
m/sec 


V 


forward velocity of rotor, m/sec 


w 

Zea. Zxo 
Zi 

4 

ttg 

Y 
6 
e 

0c 

9o 

00 

l-i 

V 

p 

a 

0 

1 

’I'o 


blade deflection at elastic axis, positive downward in 
direction normal to blade chord, m 

spanwise position of blade section from rotation axis, m 

positions of airfoil elastic and pitch axis, positive for- 
ward of midspan point, m 

quarter-chord position, positive forward of pitch axis,m 
airfoil section angle of attack, 9 + Y + 0, radians 
shaft tilt angle, positive when tilted rearward, radians 
elastic torsional rotation about elastic axis, radians 
feedback for external damping, N/m/sec 

geometric angle of attack, neglecting elastic torsional x-ota- 
tions, with respect to plane perpendicular to rotation axis, radians 

cyclic pitch amplitude, constant over span, radians 

geometric angle of attack at blade root (collective pitch), radians 

geometric linear twist, positive when 0 increases from 
root to tip, radians /m 

advance ratio, V/flR 

induced velocity, assumed constant throughout disk and 
positive downward in Yq direction (Figure 1), m/sec 

air density, kg/m^ 

solidity ratio, bc/rrR, m 

inflow angle, tan"^ ' ^^dians 

0 + Y , radians 

cyclic pitch phase angle 

angular velocity, radian/ sec 
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